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Abstract. The next-to-next-to-leading order spin-squared interaction potential for generic 
compact binaries is derived for the first time via the effective field theory for gravitating 
spinning objects in the post-Newtonian scheme. The spin-squared sector is an intricate one, 
as it requires the consideration of the point particle action beyond minimal coupling, and 
mainly involves the spin-squared worldline couplings, which are quite complex, compared 
to the worldline couplings from the minimal coupling part of the action. This sector also 
involves the linear in spin couplings, as we go up in the nonlinearity of the interaction, and 
in the loop order. Hence, there is an excessive increase in the number of Eeynman diagrams, 
of which more are higher loop ones. We provide all the Eeynman diagrams and their values. 
The beneficial “nonrelativistic gravitational” fields are employed in the computation. This 
spin-squared correction, which enters at the fourth post-Newtonian order for rapidly rotating 
compact objects, completes the conservative sector up to the fourth post-Newtonian accuracy. 
The robustness of the effective field theory for gravitating spinning objects is shown here 
once again, as demonstrated in a recent series of papers by the authors, which obtained all 
spin dependent sectors, required up to the fourth post-Newtonian accuracy. The effective 
field theory of spinning objects allows to directly obtain the equations of motion, and the 
Hamiltonians, and these will be derived for the potential obtained here in a forthcoming 
paper. 
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1 Introduction 

The anticipated direct detection of gravitational waves (GW) may become a realistic prospect 
in view of the imminent worldwide operation of second-generation ground-based interferom¬ 
eters, such as the two of Advanced LIGO [1], Advanced Virgo [2], and KAGRA [3]. This will 
transform GW science into a powerful observational tool. Inspiralling binaries of compact 
objects are considered as most promising sources for the detection of the GW signal. The 
post-Newtonian (PN) approximation of General Relativity stands out among the various and 
complementary approaches to model these systems, as it enables to an analytical treatment 
of the inspiral phase of their evolution [4]. Using the matched filtering technique to search the 
GW signal results in a pressing demand to supply accurate theoretical waveforms, where the 
continuous signal is modeled via the Effective-One-Body (EOB) approach [5]. It turns out 
that even relative high order corrections beyond Newtonian gravity, such as up to the sixth 
PN (6PN) order, are required to obtain such waveforms, and furthermore to gain information 
about the inner structure of the constituents of the binary [6]. 

Of these compact objects black holes are expected to be near-extremal, i.e. to rotate 
rapidly. Hence, considering the recent completion of the 4PN order point-mass correction 
[7], it is essential to obtain the same PN accuracy for binaries with spinning objects. To 
this end finite size effects with spins should also be taken into account. The leading order 
(LO) finite size spin effects, which are of the LO spin-squared interaction at the quadrupole 
level, and enter already at the 2PN order, were first derived for black holes in 1975 in [8], 
see also [9]. Generic quadrupoles, required to describe neutron stars, were already included 
in [9, 10], and the proportionality of the quadrupole to spin-squared was introduced in [11]. 
However, the next-to-leading order (NLO) spin-squared interaction was treated much later 
in the following series of works [12-15], and in [16] within an effective field theory (EFT) for 
spinning objects, which was formulated there. Furthermore, the LO cubic and quartic in spin 
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interactions, which enter at the 3.5PN and 4PN orders, respectively, were first approached 
in parts for black hole binaries in [14, 17]. Recently, these have been obtained for generic 
compact objects for the cubic in spin interaction in [18, 19], and also for the quartic in spin 
interaction in [18], where the latter was using the EFT for spinning objects [16]. These have 
been also partially treated in [20]. 

In this paper we derive the next-to-next-to-leading order (NNLO) spin-squared inter¬ 
action potential via the EFT for gravitating spinning objects in the PN scheme [16]. The 
NNLO level in the spin dependent case, which was already tackled using EFT techniques 
in [21], was also obtained for the spin-orbit interaction at the 3.5PN order for rapidly ro¬ 
tating compact objects [22]. The NNLO quadratic in spin potentials here and in [21] enter 
at the 4PN order. The EFT for gravitating spinning objects [16] is based on the original 
and innovative EFT approach for the binary inspiral, which was presented in [23, 24]. The 
seminal works in [25, 26] treated spin in an action approach in special and general relativity, 
respectively, and an extension of EFT techniques to the spinning case was approached in 
[27], which later adopted a Routhian approach from [28]. The EFT for gravitating spinning 
objects [16] however actually obtains an effective action at the orbital scale by integrating 
out all the field modes below this scale, and enables the relation to physical observables. 
That is, the line of work in [12, 27] leaves the temporal components of the spin tensor in the 
final results, which must be eliminated in order to obtain, e.g. the gauge invariant binding 
energy, while these components depend on field modes at the orbital scale. We show here the 
robustness of the EFT for gravitating spinning objects [16], as was already demonstrated in 
[16, 18, 22], and allows to obtain straightforwardly the physical EOM [16, 29], as well as the 
Hamiltonians from the effective action [16]. Indeed, from the potential, which we obtain here, 
we derive the EOM and corresponding NNLO spin-squared Hamiltonian in a forthcoming 
paper. 

The spin-squared sector is an intricate one, all the more so considering its LO and NLO 
levels [16], compared to the NNLO here. This sector requires the consideration of the point 
particle action beyond minimal coupling, and involves mainly contributions with the spin- 
squared worldline couplings, which are quite complex, compared to the worldline couplings, 
that follow from the minimal coupling part of the action. Yet, we also have contributions 
to this sector, involving the linear in spin couplings, as we go up in the nonlinearity of the 
interaction, and in the loop order. Thus, there is an increase in the number of Feynman 
diagrams and in complexity with respect to the NLO, which is even more excessive than 
that, which occurs in the spinl-spin2 sector [16, 21, 30]. We have here 64 diagrams, of which 
more are higher loop ones, where we have 11 two-loops, and tensor integrals up to order 
6. However, the application of the “nonrelativistic gravitational” (NRG) fields [31, 32] in 
the spinning case, as demonstrated in [16, 18, 21, 22, 30, 33], makes the computations more 
efficient. This ingredient is actually included together with other beneficial gauge choices 
in the EFT for gravitating spinning objects [16]. Finally, the EFT computations are also 
automatized here, using the xAct free packages for Mathematica [34, 35]. 

This paper is organized as follows. In section 2 we recall the required effective action 
beyond minimal coupling, following the EFT for gravitating spinning objects[16], and provide 
the Feynman rules relevant for this sector. In section 3 we provide the Feynman diagrams of 
the NNLO spin-squared interaction, and their values. In section 4 we introduce the NNLO 
spin-squared potential, and explain how to obtain from it the EOM,and the Hamiltonian. 
Finally, our conclusions are presented in 5. 
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2 The EFT for gravitating spinning objects in the PN scheme 

In this section we briefly review the effective action, which should be considered in this sector 
beyond minimal coupling, following the EFT for spinning objects, which was formulated in 
[16]. From this action the Feynman rules, which are required for the EFT computation 
of the NNLO spin-squared interaction, are derived. We build on the series of works in 
[16, 18, 21, 22, 29, 30, 33], where we follow similar notations and conventions, in particular, 
to those in [16, 18, 22]. We also use the “NRG” fields, which were already shown to perform 
well in spin dependent sectors in [16, 18, 21, 22, 30, 33]. Similarly, in what follows d is set 
to 3, as only the dependence in the generic d dimensional, related with the loop Feynman 
integrals, see appendix A in [21], should be followed. 

Let us first write the approximate metric in terms of the “NRG” fields to the orders 
required in this work: 

_ ( ^ 

_ / 1 + 2(/) -h 2(/)2 + |(^3 -Aj - 2Aj(p \ 

V ~ “^Aicj) ~^ij + — <Tij — 24 ?5ij + 2 (j)aij + AiAj + J 

( 2 . 1 ) 


The effective action for the binary system reads 

2 

5 = 5g +j;5(i)pp, (2.2) 

1=1 

where Sg is the pure gravitational action, and ^(ppp is the worldline point particle action for 
each of the two components of the binary. The gravitational action comprises the Einstein- 
Hilbert action plus the harmonic gauge-fixing term, that is 

= 5eh + 5gf = I d^x^R+ ^ J (2.3) 

where T^ = ■ 

From the gravitational action the propagators, and the self-gravitational vertices are 
derived. The “NRG” scalar, vector, and tensor held propagators, and the corresponding 
time dependent propagator vertices can be found, e.g., in section 2 of [22]. 

As for the three-graviton vertices, required for the NNLO of the spin-squared interaction, 
they are similar to those, required in [22], except for the following vertex: 

j {2aijdi(j)dj(l) - ajjdicpdicp) , (2.4) 

in which only the stationary part of the vertex is required here, unlike in the NNLO spin-orbit 
sector, where its time dependence should be included [22]. 

Also, the only four-graviton vertex, required at the NNLO level, is the stationary one, 
which is found in [21, 22]. 

We proceed to the point particle action of each of the spinning objects. The minimal 
coupling part of the covariant action, which is also invariant under rotational variables gauge 
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transformations reads [16] 



(2.5) 


with the affine parameter A, parametrized here as the coordinate time, i.e. X = t = x^, the 


4-velocity = dx^/dX, and the generic angular velocity and spin variables of the object 


respectively [16]. Notice that henceforth we drop the hat notation on the rotational 

variables. 


Yet, in the spin-squared interaction, we consider finite size effects of the objects, i.e. also 
the nonminimal coupling part of the action. Therefore, the action in eq. (2.5) should be 
supplemented here with its spin-induced part, which is quadratic in the spin [16]. From 
the symmetry considerations outlined in [16], in particular parity invariance, which holds for 
macroscopic objects in General Relativity, and the 50(3) invariance of the body-fixed spatial 
triad, it follows that the only such operator, which appears in the action [16], is given by 



( 2 . 6 ) 


where we have here the Wilson coefficient 0 ^: 52 , that is the quadrupolar deformation constant 
due to spin from [11], and similarly introduced in [12]. This coefficient equals unity in 
the black hole case. This LO nonminimal coupling for the spin-induced quadrupole, which 
also contributes the LO finite size effects, is the only operator, which contributes to the 
conservative spin-squared sector. We expect then to have in the spin-squared interaction 
both contributions from spin-squared couplings from eq. ( 2 . 6 ), and from nonlinear effects 
of the linear in spin couplings, which arise from the minimal coupling part of the action in 
eq. (2.5), that would play a smaller role in the interaction. The mass couplings also play a 
smaller role in this interaction. 

We begin then with the mass couplings, where we present their expansion in the velocity, 
to the order required here. The one-graviton couplings to the worldline mass read 



(2.7) 



( 2 . 8 ) 



(2.9) 


with similar notations as in, e.g., [ 22 ]. 

The two-graviton mass couplings, required here, read 



( 2 . 10 ) 


- 4 - 
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= m dt (j)Aiv‘^ 

1 + ... 

" ^ J 



The three-graviton mass coupling, that is require here, reads 


J dt [1 -\ -]. (2-12) 

Next, we consider the spin couplings, required for this sector, given in the canonical 
gauge [16]. Yet, we recall the kinematic terms, which involve spin with no field coupling [16], 
to the order required here 

Lkin = -S ■ O + • u X a ^1 + ) (2-13) 

with similar notations as in [22]. 

Then, the one-graviton linear in spin couplings, required in this work, read 

dt eijkSk {^diAj + [diAj - djAi) + v'-dtA^ , (2.14) 

dt 2eijkSkv"dj(j), (2.15) 

dt ^eijkSkv''diaji, (2.16) 



with similar notations as in [22]. 

The two-graviton spin coupling, which we require here, reads 



(2.17) 


Finally, we turn to the spin-squared couplings from eq. (2.6), which play the major role 
in this interaction. 

The one-graviton spin-squared couplings read 



2m 


[didjfj) ^1 -k -k - 3didk(l)v^v’^ ^1 - 

—2dtdi4>v^ + ^dkdi(j)v''v^v^v'' — dtdkcj) v'" — d‘^4>v^v^ 

- (^didifp 1^1 -k -k - didj^v^v^ ^1 - 
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= dt — 


+2dtdi4>v" ^1 + + Idfcj) (l + 

^ S^S^ (^(^d^djAkV^^ - dAAjv'^ - dtd^Aj'^ (l + 

+^ (dk {diAi - diAi) + dt {dtAk + dkAt) 

-S^ ( {dAA,v^ - - mA.) (l + IvA + 


(2.18) 


(2.19) 


= dt - 


Ces^ 

4m 


S^'S^ [didjCFkiv^v'- + dkdiCFijv'^v'- - 2didk(Tjiv^v' 


k„.l 




-2dtdiajkv^ + 2dtdk(TijV^ + d^atj^ 
{dkdkaijv'^v^ + didjakkv'v^ - 2didjaikV^v^ 
-2dtdiaijV^ + 2dtdiajjv'- + dtau)] , 


( 2 . 20 ) 


where the black square boxes represent the spin-squared operator on the worldlines, as in 
[16, 18]. 

The two-graviton spin-squared couplings, which are required here, read 






5di4> dk4> -|- 2 (5fc(/>)^ 


+2dt(l) di(j) + 3(pdidj(j) ^1 - 3(j)didkcj) + 2(f)dtdicf) 

- 5 ^ (^di(f)f ^1 -h - 3di(pdj(l)v'’v^ - 2dt(l)dicf)v'^ -4{dt(j))'^ 

+34> didicf) ^1 - -h (f)didj(l)v'‘v^ - 2(f) dtdi4> u* - 2(f)dt^ , 

( 2 . 21 ) 



Ces^ 

2m 


(^didj(f) Akv^ 


2dtdi(f) Aj 


6di(f)djAkV^ — 6di(f) dkAjV^ 


- dk(f) diAk + dk(f) dkAi - dt(f) dtAj - 4di(f) dtAj 
+3(f)didjAkV^ — 3(f>didkAjV^ — 3(f)dtdiAj'j 
+ {2dtdi(f)Ai - dkdk(f)Aiv^ + 7di(f)diAjV^ + 5di(f>djAiV^ + dt(f)diAi 

+2di(f) dtAi - 3(f) dkdkAiv’- + 3(f) didjAtv^ + 3(f) dtdiA^ , 

( 2 . 22 ) 
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= dt - 


Ces^ 

2 m 


5 * 5 ^ ( aikdjdkcp + dk4> ( diajk - ^dkcrij 


-5^ ( Uijdidjcj) + di4> ( djUij - ^dia 


jj 


= f dt ^ [S^S^ {OkAdkAj + diAkdjAk - 2d^AkdkA,) 


-2S^ [{diAjf - diAjdjAi 
Finally, the required three-graviton spin-squared coupling reads 



f 


= dt 


2 m 


(2.23) 


(2.24) 


( (I)di4)dj(l) + di dj4) \ - (cj) 


(2.25) 


Note the complexity of the spin-squared couplings with respect to the other worldline 
couplings. 


3 Next-to-next-to-leading order spin-squared interaction 

In this section we present the perturbative expansion of the NNLO spin-squared two-body 
effective action in terms of Feynman diagrams, and we provide the value for each diagram. 
At the NNLO level, which is up to order G^, the classification of the Feynman diagrams is 
similar to that in [21, 22], as can be seen in figures 1-6 below. There is a total of 64 diagrams, 
which make up the NNLO spin-squared interaction. We note that whereas the linear inter¬ 
action involves exclusively the spin-squared couplings, as we go up in the nonlinearity of the 
interaction, that is in the loop order, but even at the two-graviton exchange diagrams, we 
start to have more diagrams, which involve the linear in spin couplings. Lastly, we also have 
take into account higher order time derivatives, which appear in the 2PN order [36], and in 
up to NLO spin-orbit, and spin-squared [16] sectors. We refer the reader to [21, 22, 33] for 
further details on the generation of the Feynman diagrams from the Feynman rules, and the 
evaluation of such Feynman diagrams. In particular, we also refer to [22] for the conventions 
and notations, used here. 


3.1 One-graviton exchange 

There are 6 one-graviton exchange diagrams in this sector, which can be seen in figure I. We 
note that all diagrams here involve only the spin-squared couplings, and add up to the ones, 
which already entered in the NLO spin-squared sector [16]. Here we require tensor Fourier 
integrals of up to order 6 as in [21], see appendix A. Our treatment of time derivatives, and 
higher order time derivative terms is similar to that in [16, 22, 29]. 

These diagrams are evaluated by: 


Fig. l(al) = 


GCiEg2m2 

2mir^ 


Sf - 3(5i • nf 




GCi^g2m2 

Amir^ 


+ 2,Sivi - 2{Si • viY 
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-X- 


-XX 


(al) 


(a2) 


(a3) (bl) 


--X-- 


(b2) 


(c) 


Figure 1. NNLO spin-squared one-graviton exchange Feynman diagrams. 


-|- 6Si ■ nvi ■ nSi ■ vi — • n)^ — 9^2 (5i • n)^ — 6Sf{vi ■ n)^ 

25i • Sivi ■ fi — Si ■ nSi ■ vi 


I GCiEs^rn2 ^ ^ ^ GCiEs^rn2 

H- 7 ,—oi • aiiii • n — 


mir^ 

— Si ■ nSi ■ vi 


mir^ 

_l_ GCiES'^'m2 ^2 _|_ GGiEs^rn2 

mir ^ IGmir^ 


SOStvfvi - 12vi{Si-viY 


— 12v2{Si ■ viY + l9Sfvf + 7Sfv2 + 4,2Si ■ nvi ■ nSi ■ vivf 

+ 36^1 • nvi ■ nSi ■ V 1 V 2 — bivfv^iSi ■ n)^ — 6('Ui • n)^(5i • viY 

— 36Sfvl{vi ■ fiY — 36Sfv2{vi ■ n)^ — 21(5i • n^vf — 21{Si ■ 


+ 


GGieS^1^2 


2mir'^ 

+ 3^1 • n^i • diV2 


Slv^di ■ n — 2vi ■ fiSi ■ viSi ■ Si -|- 2Sfvi ■ nvi ■ Si 
GGies^1^2 


2vi ■ nSi ■ viSi ■ vi — 2Si ■ SiVi ■ nvi 


2mir‘^ 

+ %Si ■ SiVi ■ nv2 — 3Si ■ nSi ■ viv\ — 3Si ■ nSi ■ V 1 V 2 
GGiES‘2'm'2 


+ 


2mir 


2Si • • Si + (2Si • viSi ■ VI - 2Si ■ Sivf + vfSf + 3vlSl) 


+ 2{Si ■ SiY + 4(5i • viSi ■ Si + Si ■ viSi ■ Si — 2Si ■ SiVi ■ Si+ vi ■ SiS^) 


Fig. I(a2) = 


- 2{Sfvl - (Si ■ viY) 

GCiES2n^2 
Amir^ 


(3.1) 


25i • FiS*! • V 2 — Sivi ■ V 2 + 3Sfvi ■ nv 2 ■ n — 65i • nSi ■ V 1 V 2 ■ n 


— OiSi • nvi ■ nSi ■ V 2 — 3Fi • F 2 (<S'i • n)^ -I- 15Fi • nv 2 ■ n{Si ■ Sy 


+ 


GCiEs^'m‘2 

2mir'^ 

GGies‘2'^2 

8mir^ 


Si ■ S 1 V 2 ■ n — Si ■ nSi ■ V 2 — Si ■ nSi ■ V 2 + 3S'i • nSi ■ nv 2 ■ n 
45i • viv^Si ■ V 2 + 3Sfv‘fvi ■ V 2 + 65i • Fi5i • V 2 V 2 — 3Sfvi ■ V 2 V 2 


— 4Fi • F2 (.Si • Fi)^ — 95^Fi • nv‘{v 2 ■ n — 12Si ■ nSi ■ viv‘lv 2 ■ n 

— 185i • nFi • nv\Si ■ V 2 + 125i • nFi • nSi ■ vivi ■ V 2 + 95fFi • nv 2 ■ nv 2 

— 185i • nSi ■ V 1 V 2 ■ nv2 — I8.S1 • nFi • nSi ■ V 2 V 2 — 9viVi ■ V2{Si ■ n)^ 

— 9Fi • F2f|(5i • sy + 12Fi • Sv2 ■ S{Si ■ Fi)^ -|- 65i • Fi5i • F2(Fi • n)^ 

— 185fFi • F2 (Fi • sy + 45Fi • Sv‘lv2 ■ S{Si ■ n)^ -|- 45Fi • Sv2 ■ Sv2{Si ■ Sy 


- 8 - 




















































+ 30Sfv2 ■ n{vi ■ n)^ — SOS'! • nSi ■ viV2 ■ n{vi ■ nY 


GCi^g2m2 


Sfvi ■ diV2 ■ n + Si ■ vidi ■ fiSi ■ V2 — vi ■ fiSi ■ diSi ■ V2 


4mir^ 

— 6Si ■ nvi ■ aiSi ■ V2 — 2Sfai ■ nvi ■ V2 — Si ■ nSi ■ aivi ■ V2 — 2Sfvi ■ nai ■ V2 
+ Si • nSi ■ viSi ■ V2 — SSfvi ■ nv2 ■ 02 + 6 S 1 • nSi • V 1 V 2 ■ 0.2 
+ GSfvi ■ ndi ■ nv 2 ■ n — 3Si • nSi ■ viSi ■ nv 2 ■ n + 35i • nvi ■ nSi ■ aiV 2 ■ n 
+ 9tTi • aii ;2 • n{Si ■ nY — Qui • nv 2 ■ 02(^1 • nY 


+ 


GCieS^1^2 

Amir'^ 


2Si ■ viSi ■ V 1 V 2 ■ n — 55i • Siv\v 2 ■ n + vi ■ nSi ■ viSi ■ V 2 


+ 3Si • nvfSi ■ V2 + vi ■ nSi ■ viSi ■ V 2 + 3Si ■ nviSi ■ V 2 

— 4Si • SiVi ■ nvi ■ V2 + Si ■ nSi ■ vivi ■ V 2 + Si ■ nSi ■ vivi ■ V 2 — 3Si • 51^2 • nv2 
+ 3Si • nSi • V 2 V 2 + 35i • nSi • V 2 V 2 — 3Si ■ nvi ■ nSi ■ V 1 V 2 ■ n 

— 3Si ■ nvi ■ nSi ■ V 1 V 2 ■ n — 9Si • nSi • nviV 2 ■ n — 9Si • nSi • nv 2 ■ 

+ 6 S 1 • S 1 V 2 ■ n{vi ■ nY - —~ ■ ^ 1*^1 • '^2 — Si • nSi • aiV 2 ■ n) 


2 mir 

+ (Si • -iTiSi • -i;2 + Si • -iJiSi • V2 - 2Si ■ SiVi ■V2-vi- V 2 SI + 2Si ■ SiVi ■ nv2 ■ ft 
— Si ■ nSi ■ V 1 V 2 ■ n — Si ■ nSi • V 1 V 2 ■ n + vi ■ nv 2 ■ nSf) + (2Si • aiSi • V 2 
+ 2Si • aiSi • V2 - 4Si • Sioi • V2 + 3Si • Si'U2 • ^2 + oi • V 2 SI 
+ 4Si • Siai • nv 2 ■ n — 2Si ■ nSi ■ diV 2 ■ n — 2Si ■ nSi ■ diV 2 ■ n 
+ 3Si • nSi • nv 2 ■ 0.2 — di ■ nv 2 ■ nSf) + 2(Si • viSi ■ V 2 — Sfvi ■ V 2 

+ Sfvi ■ nv2 • n — Si • nSi • 'Uii;2 • n) 


, GCiEs2m2 ® - - 

H-^-Si t>2 • n 


2 mi 


(3.2) 


Fig. I(a3) = 


8S1 • FiSi • V2V1 ■ V 2 - 3 S(viv 2 + 2 i; 2 (Si • viY + 2 vi{Si ■ 7/2)^ 


GGiES‘2ni2 
16mir^ 

— 6Si(i;i • V 2 Y — 2Avi ■ nSi ■ V 1 V 2 ■ nSi ■ V 2 — 12Si • nv‘fv2 ■ nSi ■ V 2 

+ 36S^'i7i • nv2 ■ nvi ■ V 2 — 24Si • nSi • V 1 V 2 ■ nvi ■ V 2 — 24Si • nvi ■ nSi ■ V 2 V 1 ■ V 2 

— 12Si • nvi ■ nSi ■ viv^ — 3vlv2{Si ■ nY + 9Sfv2{vi ■ nY — 6(Si • viY{v2 ■ nY 
+ 9Sjvl{v2 ■ nY - 6(Fi • re)^(Si • V 2 Y - 6(Si • nY{vi ■ V 2 Y 

+ GOFi • nv2 ■ nvi ■ V2{Si ■ nY + 6OS1 • nv2 ■ nSi ■ V2{vi ■ nY 
+ 15^2(Si • nY{vi ■ nY — 45Sf(Fi • nY{v2 ■ nY + 6OS1 • nvi ■ nSi ■ vi{v2 ■ nY 
+ 15ff(Si • nY{v2 • nY — 105(Si • nY{vi ■ nY{v2 ■ nY 


+ 


GGiEs^ni2 

IGmir^ 


4 Si • aii;2 • nSi ■ V2 — 6 Sfv 2 ■ ndi ■ V2 + 4 Si • nSi ■ V 2 di ■ V2 


— 3S‘ldi ■ nv 2 + 2Si ■ nSi ■ aiff + 3Sifi02 • n — 4'Ui • nSi ■ viSi ■ 02 

— 2 Si • fivlSi ■ (12 + GS^Fi • fivi ■ 02 — 4 Si • nSi • viVi ■ 02 — 2 d 2 ■ n(Si • viY 
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Fig. l(bl) 


+ 2 ai ■ n{Si ■ V 2 )^ — 125i • nai ■ nv 2 ■ nSi ■ V 2 + 125i • nvi ■ nSi ■ via 2 ■ n 
— 6 v 2 ■ nai ■ V2{Si ■ n)^ — 3oi • ■ n)^ + 3v‘fa2 ■ n{Si ■ re)^ 

+ 6 Fi • nvi ■ a2{Si ■ n)^ — 9Sfa2 ■ n{vi ■ n)^ + 65i • nSi ■ a2{vi ■ n)^ 

+ DiS^ai • n{v2 ■ n)^ — 65i • n5i • di{v2 ■ n)"^ — 15a2 • n{Si ■ n)‘^{vi ■ n)^ 


+ 15ai • n{Si ■ n)‘^{v2 ■ nY 


+ 


GCi^g2m2 


Amir‘s 

+ 251 • V 1 V 2 ■ nSi ■ V2 + 2 Fi • n5i • ?;25i • V2 


2Si ■ V 1 V 2 ■ fiSi ■ V2 
— 65i • 5 i?;2 • nvi ■ V 2 


+ 25i • n5i • V 2 V 1 ■ V2 + 2Si ■ nSi ■ V 2 V 1 ■ V2 — 35i • + 5i • n5i • V 1 V 2 


+ 5i • nSi ■ V 1 V 2 — 65i • nFi • nv2 ■ nSi ■ V 2 — 65i • nFi • nv2 ■ nSi ■ V 2 

— 65i • n5i • nv2 ■ nvi ■ V 2 — 35i • n5i • nvi ■ nv^ + 95i • SiVi ■ n{v2 ■ nY 

— 3Si ■ nSi ■ vi{v2 ■ nY — 3Si ■ nSi ■ vi{v2 ■ nY + 155i • n5i • nvi ■ n{v2 ■ nY 


GCiEg2m2 

16mir 


2(25i • 'i;25i • V2 — 35i • S 1 V 2 — 2Si ■ nv2 ■ nSi ■ V 2 


— 2Si ■ nv2 ■ nSi • F 2 — 5i • nSi ■ nv2 + 35i • 5i(F2 • nY + 35i • nSi ■ n{v2 ■ nY) 

— ( 2 S 1 ■ aiSi ■ 0,2 — 35iai • 0,2 + 35iai • na 2 ■ n — 2Si ■ nSi ■ 010,2 ■ n 

— 2 Si ■ noi ■ nSi ■ 02 — oi ■ 02(81 ■ nY + 3 ai • na2 ■ n(Si ■ nY) — 4 ( 5 i • • 02 

+ 5 i • viSi ■ 02 — 35 i • 5 iFi • 0,2 + 35 i • SiVi ■ nd2 ■ ft — Si ■ nSi ■ vid2 ■ n 

— Si ■ nSi ■ vid2 • n — 5 i • nvi ■ nSi ■ 0,2 — Si ■ nvi ■ nSi ■ 02 — Si ■ nSi ■ nvi ■ 02 


+ 35 i • n 5 i • nvi ■ 00,2 ■ ft) — 2 ( 35 iu| — 2 ( 5 i • V2Y + 45 i • nv2 ■ nSi ■ V2 
+ vl(Si ■ nY - 3 SI(v 2 ■ nY - 3(5i • nY(v2 ■ nY) 


GCiEg2m2 

8nii 


(35i • 5 ia 2 • n — 5i • n5i ■ 02 — Si ■ nSi ■ 02 + Si ■ nSi ■ na 2 ■ n) 


+ (35i02 • n — 2Si ■ nSi ■ 02 + 02 ■ n(Si ■ nY) 

2GCiES^ni2 \ „2^ Q- - /c -^2 , 

b^vi ■ V 2 - 3vi ■ V 2 {Si ■ n) + 

GCiES‘^ni2 


(3.3) 


mir^ 

— Si ■ nSi ■ V2 — Si ■ nSi ■ V2 


2GCiEs^ni2 


2Si ■ S 1 V 2 ■ n 


mir-^ 


mir^ 

o o2 2 —* —* I 02 —* —* 2 
3Dit;iUi • t ;2 + 5i?;i • V 2 V 2 


— 2vi ■ V2(Si ■ viY + 65i • nFi • n5i • vivi ■ V2 — 3v\vi ■ V2{Si ■ fty 


— 3vi ■ V2V2 {Si ■ nY — QSfvi ■ V2(vi ■ fty 


+ 


GGiEs'^nfi2 


niir^ 


2S‘lvi ■ diV 2 ■ n 


+ 5i • Fiai • n5i • 'i;2 + Fi • nSi ■ oiSi ■ V 2 — 2Si ■ nvi ■ oiSi ■ V 2 

— 25 iai • nvi ■ V2 + Si ■ nSi ■ oivi ■ V2 — 2 S‘fvi ■ noi ■ V2 + Si ■ nSi ■ viOi ■ V2 


+ 


GGiEs^ni2 


mir^ 


2Si ■ Sivlv 2 ■ n + vi ■ nSi ■ Fi5i • F 2 — 5i • nvlSi ■ V 2 


+ Fi • nSi ■ Fi5i • F2 — 5i • nvISi ■ V2 — 45i • 5iFi • nFi ■ V2 + Si ■ nSi ■ viVi ■ V2 
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+ Si • nSi ■ vivi ■ V 2 + 25i • S1V2 ■ n?;| — Si • nSi • V2V2 — Si • nSi • V2V2 (3.4) 


Fig. I(b2) = - 


GCi^g2m2 


mir'^ 


2Si • FiSi • V2V1 ■ V 2 — Si{vi ■ V 2 )^ + SSfFi • nv 2 ■ nvi ■ V 2 


— 6S1 • nSi • V1V2 ■ nvi ■ V 2 — 6S1 • nvi ■ nSi • V2V1 ■ V2 — 3(Si • n)‘^{vi ■ V 2 y 


+ 15i;i • nv2 ■ nvi ■ V2{Si ■ nY 


+ 


GCiEs'^'fn2 

mir^ 


81^2 ■ nai ■ V2 


— 2 Si • nSi • i 72 ai • V2 — Sfvi ■ nvi ■ 0,2 + 2 Si • nSi • v\V\ ■ 02 
+ 3F2 • nai ■ V2{Si ■ nY — 3 'Ui • nvi ■ 02(Si • nY 


GCies‘^^2 


mir^ 


Si • 'Ui'U 2 • nSi ■ V2 + Si ■ V 1 V 2 ■ nSi ■ V 2 + 2vi ■ nSi ■ V 2 S 1 ■ V 2 


— 6S1 • S1V2 ■ nvi ■ V2 + 3 Si • nSi • i;2i;i • V2 + 3 Si • nSi • V2V1 ■ V2 

— 2 Si • SiFi • nv2 — 3 Si • nvi ■ nv2 ■ nSi ■ V2 — 3 Si • nvi ■ nv2 ■ nSi ■ V2 

— 6S1 • nSi • nv2 ■ nvi ■ V2 + 6S1 • SiFi • n{v2 ■ nY 
GGiEs^n%2 


+ 


mir 


(2S1 • F2S1 • V 2 — 2Si • S1U2 — Si • nv 2 ■ nSi ■ V2 


— Si • nv2 ■ nSi ■ V2 + 2 Si • Si(F 2 • nY) + (Siai • 02 + oi • a 2 (Si • n)^) 

— (Si • FiSi • a2 + Si • FiSi • 0,2 — 4 Si • SiFi • 0,2 + 2 Si • SiFi • 00,2 ■ n 

— Si • fivi ■ nSi ■ 02 — Si ■ nvi ■ nSi ■ 02 — 2Si • nSi • nvi ■ 02 ) — 2[SiV2 


- (Si • V 2 )'^ + Si • nv 2 ■ nSi ■ V2 - Sf (F2 • nY) 


+ 


GCiEs'2ni2 


nil 


Fig. 1(c) =- 


— Si • nSi • a 2 — Si • nSi • 02 ) + 2 (Sia 2 ■ n — Si ■ nSi ■ 02 ) 
GCiES‘2m2 


( 2 S 1 ■ §102 ■ n 

(3.5) 


miv^ 

GCiES‘^ni2 

niir"^ 


S1V1V2 - Si{vi ■ V2) - 3 viV 2 {Si-n) + 3(Si • re) (Fi • F2)" 
2S1F2 • real • V2 — 2Si • reSi • ^201 • V2 — Sfoi ■ reu^ 


+ 2Si • reSi • 01^2 + ai • re(Si • ^2)^ 


2GCiES‘^nri2 


mir^ 


2 Si ■ S1V2 ■ nvi ■ V 2 


— Si ■ reSi • V2V1 ■ V2 — Si ■ reSi • V2V1 ■ V2 — 2 Si ■ SiVi ■ rev^ + Si • reSi • Fiu| 


+ Si • reSi • V1V2 


2GCiEs^ni2 


mir 


(Si • V2S1 • F2 + Si • Sivl) 


+ ('^ 1^2 4" ('S'l ■ ^ 2 )^) 


(3.6) 


3.2 Two-graviton exchange and cubic self-interaction 
3.2.1 Two-graviton exchange 

There are 12 two-graviton exchange diagrams, which can be seen in figure 2. They involve 
two-graviton mass couplings, which at the NNLO level are also contracted with linear in spin 
couplings, or two-graviton spin-squared couplings. We note that time derivatives originate 
here also from the spin-squared couplings. 
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Figure 2. NNLO spin-squared two-graviton exchange Feynman diagrams. 


These 
Fig. 2(al) 


diagrams are evaluated by: 
2 G^m2 


Sfvf — {Si ■ viY + 2Si ■ nvi ■ nSi ■ vi — vi{Si ■ ny — Syvi ■ fty 


(3.7) 


Fig. 2(a2) = - 


4G^m2 


5i • viSi ■ V 2 — Sivi ■ V 2 + Sivi ■ nv 2 • n — 5i • nSi ■ viV 2 ■ n 


— Si ■ nvi ■ nSi ■ V2 + vi ■ V2{Si ■ fty 


Fig. 2(bl) = - 


G‘^Cies^^2 


Sf - 3(5i • nf 


G‘^Cies^'^2 

4j,4 


(3.8) 

8Sjvf - 9Sjv^ - 2{Si ■ viY 


+ 6Si ■ nvi ■ nSi ■ vi — 18vf{Si ■ fiy -|- 27t;|(5i • n)^ — 6Sl{vi ■ fty 

2Si ■ SiVi ■ n — Si ■ nSi ■ vi 


G‘^GiEs^'m‘2 ^ ^ ^ ^ G'^C'iES2m2 

-^-Di • oioi • n -|-^- 


Fig. 2(b2) = - 


— Si ■ nSi ■ vi 

G‘^Cies^'<^2 


G'^Gies^'itT'2 y 


Sj 


(3.9) 


4^4 


2S'i • FiS*! • V 2 — SiVi ■ V 2 + 25i • nvi ■ nSi ■ vi -|- 4S'^Fi • nv 2 ■ n 


— 8Si ■ nSi ■ viV2 ■ ft — 65i • nvi ■ nSi ■ V2 — 3Fi • F2(5i • n)^ — 5'^(Fi • nY 


+ 18Fi • nv2 ■ n{Si ■ fiy — 3(S'i • n)^(Fi • fty 


G'^Gies'^1T^2 


Si ■ SiVi ■ n 


— 2 Si ■ SiV2 ■ n + Si ■ nSi ■ V2 + Si ■ nSi ■ V2 + Si ■ nSi ■ nvi ■ n 

— 45 i • nSi ■ nv2 ■ n 


(3.10) 


Fig. 2(b3) = - 


G'^Gies^'1^2 

47-4 


Sfvi ■ V2 + OiSi • fivi ■ fiSi ■ vi — ^S\vi ■ nv2 ■ ft 


— 65i • nvi ■ nSi ■ V 2 — 3Fi • F2(5i • n)^ -|- 3iS'7(Fi • n)^ -|- 18Fi • nv 2 ■ n{Si ■ nY 


— 15(5i • ny{vi ■ fiy 


+ 


G‘^Gies‘^'^2 

27-3 


Si ■ iSiFi • fi — 35i • fiSi ■ fivi ■ fi 


(3.11) 
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Fig. 2(b4) 
Fig. 2(b5) 


Fig. 2(cl) 


Fig. 2(c2) 


Fig. 2(c3) 


Fig. 2(c4) 


Fig. 2(c5) 


2G*^C*iEg2W?-2 

2 G'^Cies‘^i ^2 


Sfvi ■ V2 - 3Fi • V2{Si ■ n)^ 
Sfvi ■ V2 - 3vi ■ V2{Si ■ n)^ 


(3.12) 


2 G'^CiEg 2 m 2 


2 Si ■ S 1 V 2 ■ n 


— Si ■ nSi ■ V 2 — Si ■ nSi ■ V 2 


(3.13) 


2 G'^Cies2'<^2 


mir^ 


Sl - 3{Si • nf 


+ 


G^G 




2mi 


Slvl + 2Si ■ viSi ■ V 2 


+ 25fFi • V 2 — 14:SiV 2 — {Si ■ Fi)^ + 45i • nvi ■ nSi ■ vi — 8SiVi ■ nv 2 ■ n 
— 8Si ■ nSi ■ V 1 V 2 ■ fi — 2vl{Si ■ n)^ + 36f|(5i • n)^ + 10Sf{v2 ■ n)" 

, G'^CiES^'ml^ -.q 2 

H-o— - 0-2 ■ nSl (3-14) 


mir 
2 mir‘^ 


6 Si ■ viSi ■ V 2 + Sfv 2 — 3(5i • ^ 2 )^ — 245i • nSi ■ viV 2 ■ n 


— 24Si ■ nvi ■ nSi ■ V 2 + 24:Si ■ nv 2 ■ nSi ■ V 2 — 12Fi • F2(5i • n)^ + 2S‘f{v2 ■ n)^ 

+ 72Fi • nv 2 ■ n{Si ■ n)^ — 36(5i • n)'^{v 2 ■ F)^ — —- Si ■ nSi ■ V 2 

J mir^ L 

+ Si ■ fiSi ■ V 2 — 4:Si ■ nSi ■ nv 2 ■ n (3.15) 

G‘^GiES^fn 2 " 


mir 


5Sfv2 — 3(5i • ^2)^ + 85i • nv2 ■ nSi ■ V 2 — 8^2 (<Si • ny 


- 8 SI{v 2 ■ ny 

2 G‘^GiEs‘^nfi 
mir^ 


(3.16) 


2 , 


45i • viSi ■ V 2 — 12S‘fvi ■ V 2 + 5 Sfv 2 — 5{Si ■ F 2 )" 


+ lASlvi ■ nv 2 ■ n — Si ■ nSi ■ V 1 V 2 ■ n — 155i • nvi ■ nSi ■ V 2 
+ 205i • nv 2 • nSi ■ V 2 + I 8 F 1 • F 2 (.Si • n)^ — 18£'^(F2 • n)^ 
2 G‘^GiEs^nT,l\ 2 ^ ^ -c - 

H- 5 - 5 D]^a 2 • n — 7oi • nbi ■ 02 


mir-^ 
2 G'^CiEs^ni 2 " 


mir^ 

- 2 SI{v 2 ■ n)^ 


2Sfv2 — {Si ■ ^2)^ + 25i • nv2 ■ nSi ■ V2 — V2{^i ' n)^ 


(3.17) 


(3.18) 


3.2.2 Cubic self-interaction 

There are 27 cubic self-interaction diagrams in this sector, which can be seen in figure 3. 
Tensor one-loop integrals up to order 4, which appear in appendix A of [21], are required to 
be initially applied here in addition. 

These diagrams are evaluated by: 


T-- G‘^m 2 ,^ ^,2 , G^m 2 

Fig. 3(al) =—^{Si • n) + 


iSfvi — 4(5i • Fi)^ -|- 105i • nvi ■ nSi ■ Fi — 2 i;f(5i • fty 


-|- 3 ^ 2(51 • fiy — 85^ (Fi • fty 




2 G‘^m 2 


Slai ■ n — Si ■ fiSi ■ ai 
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(dl) 


Fig. 3(a2) 


Fig. 3(a3) 



Figure 3. NNLO spin-squared one-loop Feynman diagrams. 


+ 


2G‘^m2 

j-Z 

G‘^m2 


Si ■ SiVi ■ n — Si ■ fiSi ■ vi 


SSfvi — 65i • viSi ■ V2 — 2SiVi ■ V2 — 12S'i • nvi ■ nSi ■ vi 


(3.19) 


-|- SSfvi ■ nv 2 ■ n + 245i • nSi ■ viV 2 ■ n + IGSi • nvi ■ nSi ■ V 2 
— 137;f (Si • n)^ -|- 8Fi • ^ 2(^1 • n)^ — 55i (Fi • n)^ — 48Fi • nv2 ■ n{Si ■ fiy 


+ 27(Si • n)^(Fi • fty 


G^m 2 

2r3 


6S1 • SiFi • n -|- 25i • nSi ■ vi 


+ 8Si ■ nSi ■ vi — 6S1 • 5 iF 2 • n — 5Si ■ nSi ■ V 2 — 55i • nSi • V 2 

G‘^m 2 


— 205i • nSi • nvi ■ n -|- 20Si • nSi ■ nv 2 ■ ft 


+ 


2 r2 


3Sf - 5(Si • fiy (3.20) 


G^m2 

2 r^ 


Si ■ viSi ■ V 2 + Sfvi ■ V 2 — 4SiFi • nv 2 ■ n — 4Si • nSi • F1F2 • n 


— 4Si • nvi ■ nSi ■ V 2 — 2vi ■ V 2 {Si ■ fty -|- 12Fi • nv 2 ■ n{Si ■ fty 
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Fig. 3(a4) 
Fig. 3(a5) 
Fig. 3(a6) 


Fig. 3(a7) 


Fig. 3(a8) 


Fig. 3(a9) 

Fig. 3(bl) 

Fig. 3(b2) 

Fig. 3(b3) 

Fig. 3(b4) 

Fig. 3(b5) 


G^m 2 

2^3 


25i • S 1 V 2 ■ n + Si ■ nSi ■ V 2 + Si ■ nSi ■ V 2 — • nSi ■ nv 2 ■ n 


4G^m2 


Si ■ viSi ■ V 2 — Sfvi ■ V 2 + 2Sfvi ■ nv 2 ■ n — 2Si ■ nSi ■ V 1 V 2 ■ n 


8G^C'iE52m2 


Slvl - 2 vl{Si ■ n)‘ 


m\r^ 
4G^C'iE52m2 

- 2 >vi ■ V 2 {Si ■ fty 


(3.21) 

(3.22) 

(3.23) 


S'lvi ■ V 2 — 2 Si ■ nSi ■ V 1 V 2 ■ n + 25i • nvi ■ nSi ■ V 2 


16G^C 


lES 2 m 2 


2 r 


mir* 

— Si ■ nSi ■ V 2 — Si ■ fiSi ■ V 2 
AG‘^GieS'2'IT^2 


sfvi ■ V2 - 2 vi ■ V2{Si ■ fty 


+ 


8G^Cigg2m2 ■■ 

mir^ 


(3.24) 
2^1 • 5ii;2 • n 

(3.25) 


'S'lFi • V 2 + 2Si ■ fiSi ■ V 1 V 2 ■ n — 2Si ■ nvi ■ nSi ■ V 2 


- 3vi ■ V 2 iSi ■ n)^ 

— 35i • nSi ■ V2 
AG‘^CiEs‘2m2 


+ 


2 G‘^Gies 2^2 


2Si ■ S 1 V 2 • n — 3Si ■ fiSi ■ V 2 


Sfvj - 3vj{Si • ny 

— 3Si ■ fiSi ■ vi — 3Si ■ fiSi ■ vi 
G‘^m 2 


2G‘^CiEs^nT’2 


2Si ■ SiVi ■ n 


3Sfvi — 3{Si ■ vi)‘^ + 105i • nvi ■ nSi ■ vi — 2vf{Si ■ fty 


— 8 Sl{vi ■ n)^ 
4G'^CiES2m2 


mir^ 


{Si ■ ^ 2 )^ — 8 Si ■ nv 2 ■ nSi ■ V 2 — 2^2(<Si • ny 


+ 12(5i • ri)^(F 2 • re)" 


G‘^GiEs^ni 2 


Sfvf + 2(5i • Fi)^ — 185i • nvi ■ nSi ■ vi — 5vi{Si ■ fty 


— 4:Sl{vi ■ re)^ + 30(^1 • re)^(Fi • re)" 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


(3.30) 


G‘^GiEs^ni 2 

4^4 


Sfv^ + 2(5i • ^ 2 )^ — 125i • nv 2 ■ nSi ■ V 2 — 5f|(Si • re)" 


- Sl{v 2 ■ fty^ + 15(5i • n)‘^{v 2 ■ fiy^ (3.31) 

+ (*^1 ■ ^ 1 )^ “ “ ^v\{Si ■ n)^ 


G‘^CiEs^ni 2 


2 mir^ 

+ 12(5i • re)^(Fi • re)^ 
G‘^CiEs'^ni 2 




5i • fiSi • oi — oi • n{Si ■ re)" 


+ 


2 mir^ 


2Si ■ SiVi ■ n — 3iS'i • nSi ■ vi — 35'i • nSi ■ vi 
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Fig. 3(b6) 


Fig. 3(cl) 


Fig. 3(c2) 


Fig. 3(c3) 


Fig. 3(c4) 


Fig. 3(c5) 


+ 85i • fiSi ■ nvi ■ n 

+ (5? + (^1 • nf) 
G‘^CiEs‘^m 2 




2 r 


2 mir^ 


{Si -31 +Si- nSi ■ n) 


(3.32) 


{Si ■ vi)"^ — 6 Si ■ nvi ■ fiSi ■ vi — vf{Si ■ n)^ + 6(5i • n)^{vi ■ fty 


G ‘^ Gies ‘^'^2 


G'^GiEs' 2 m 2 


2 r^ 

+ 165i • n5i • nvi ■ ft 
G‘^m 2 


Slai ■ n — Si ■ nSi ■ oi + 2oi • n{Si ■ fty 


6 Si ■ SiVi ■ n — 5Si ■ nSi ■ vi — 5Si ■ nSi ■ vi 


+ 


AG'^Gies‘^'1^2 


Si ■ nSi ■ n + {Si ■ fiy 
2S\vl — Si ■ viSi ■ V2 + S\vi ■ V2 — 2(5i • Fi)^ + 65i • nvi ■ nSi ■ vi 


(3.33) 


— ^S\vi ■ nv 2 ■ n + 45i • ftSi ■ viV 2 ■ n — 2 vl{Si ■ fty — 4:Sf{vi ■ fty 


+ 


G‘^m 2 


Sfdi ■ fi — Si ■ nSi ■ ai 


+ 


G‘^m 2 ' 


Si ■ SiVi ■ fi — Si ■ nSi ■ vi (3.34) 


AG^GiEs^mi 


2 , 


mir^ 


Si ■ viSi ■ V 2 + {Si ■ V 2 y - 4^1 • nSi ■ viV 2 ■ ft 


— 45i • nvi ■ nSi ■ V 2 — 8 Si ■ nv 2 ■ nSi ■ V 2 — 2vi ■ V 2 {Si ■ fiy 

— 2t>|(5i • fiy + 12Fi • nv 2 ■ fi{Si ■ fiy + 12(5i • n)^(F2 • fiy 

Si ■ fiSi ■ V 2 + Si ■ fiSi ■ V 2 — 4Si • fiSi ■ fiv 2 ■ fi 


AG'^GiEs^'iTy ^' 


mir- 


(3.35) 


G‘^Gies ^'^2 


Slvl + 2Si ■ viSi ■ V 2 + Slvi ■ V 2 + 2{Si ■ Fi)' 


— 185i • nFi • n5i • Fi — 45fFi • fiv 2 ■ n — 8Si ■ nSi ■ F 1 F 2 • n 

— 105i • nFi • fiSi ■ V 2 — 5fi(5i • n)^ — 5Fi • F 2 (.Si • n)^ — 45i(Fi • n)^ 

G‘^Cies' 2 ' 1^2 


+ 30Fi • fiv 2 ■ fi{Si ■ fiy + 30(5i • F)^(Fi • fiy 


+ 35i • fiSi ■ Fi + 35i • fiSi ■ Fi — 125i • nSi ■ nvi ■ n 


2Si ■ SiVi ■ fi 

(3.36) 


G ‘^ Gies ‘^'^2 

" 2^^4 


25i • Fi5i • V 2 + SiVi ■ V 2 — SiVi ■ fiv 2 ■ n — OSi • nSi ■ F 1 F 2 • n 


— 65i • FFi • fiSi ■ V 2 — 5Fi • F2(5i • n)^ + 15Fi • fiv 2 ■ fi{Si ■ fiy 


+ 


G'^GiEs^nT'2 


Si ■ fiSi • F 2 + 5i • fiSi ■ V 2 — 2 Si ■ fiSi ■ fiv 2 ■ fi 


(3.37) 


G GiES 2 m 2 „2- ^ AC ^ ^ ^ 

- -^• viSi ■ V2 + SiVi ■ V2 - 4 ii • nSi ■ V1V2 ■ n 

mir ^ L 

— 45i • FFi • fiSi ■ V 2 — 4Fi • F2(5i • n)^ + 12Fi • nv 2 ■ fi{Si ■ n)"^ 
2Si ■ S 1 V 2 ■ fi — 3Si ■ fiSi ■ V 2 — 8 Si ■ fiSi ■ V 2 


G'^GiEs’^ny 


2 mir^ 


- 16 


















































(3.38) 


Fig. 3(c6) 

Fig. 3(dl) 

Fig. 3(d2) 

Fig. 3(e) 

Fig. 3(f) 

Fig. 3(gl) 

Fig. 3(g2) 


+ 85i • fiSi ■ nv 2 ■ n 

- ^ 

r* L 

— 45i • nSi ■ viV 2 ■ n — 2 Si ■ nvi ■ nSi ■ V 2 — vfiSi ■ n)"^ — vi ■ V 2 {Si ■ fty 


+ 6Fi • nv 2 ■ n{Si ■ n)^ + 6(5i • n)^(Fi • fty 


+ 


G‘^CiEs‘2'm2 

2^3 


2 Si ■ Sivi ■ n 


— 35i • nSi ■ vi — 35i • nSi ■ vi + 45i • £'iF 2 • n — 2 Si ■ nSi ■ V 2 

— 2Si ■ nSi ■ V 2 + 8 Si ■ nSi ■ nvi ■ n + 85i • nSi ■ nv 2 ■ n 


G^m2 


Sfv2 — 45i • nv2 ■ fiSi ■ V2 — 3?;2(*S'i • n)^ — 35jf(F2 • n)" 


4^4 

+ 9{Si ■ ny{v2 ■ fty 

G'^m 2 


(3.39) 


(3.40) 


35i • Fi5i • V 2 — 35'^Fi • V 2 + lOS’fFi • nv 2 ■ n — 65i • nSi ■ viV 2 ■ n 


— ASi ■ nvi ■ nSi ■ V2 + 2Fi • F2(5i • n)^ 

G‘^m 2 


(3.41) 


SSivf + 2Si ■ viSi ■ V 2 — 2Sivi ■ V 2 — 3(S'i • Fi)^ + 145i • nvi ■ nSi ■ vi 
+ SSfvi ■ nv 2 ■ n — 8 Si ■ nSi ■ viV 2 ■ n — 45i • nvi ■ ftSi ■ V 2 
— 2 vi ■ V 2 {Si ■ fty — 8S'^(Fi • rVy + 12Fi • nv 2 ■ n{Si ■ fty — 12(5i • n)^(Fi • n)^ 


+ 


G‘^m 2 


35i • SiVi ■ fi — 2Si ■ nSi ■ vi — Si ■ nSi ■ vi + 45i • nSi ■ nvi ■ n 


(3.42) 


G‘^m 2 

2r4 


4£'i • viSi ■ V 2 — 5Sfvi ■ V 2 — ‘ 2 S 1 V 2 + 2(S'i • ^ 2 )^ + 7S'fFi • nv 2 ■ ft 


— 2Si ■ nSi ■ V 1 V 2 ■ fi — 65i • fivi ■ nSi ■ V 2 — 85i • fiv 2 ■ nSi ■ V 2 


+ 7vi ■ V 2 {Si ■ fiy + 85'^(F 2 • fty — 9Fi • nv 2 ■ n{Si ■ fty 


+ 


G‘^m 2 


£'^02 • n 


2 G‘^m 2 y ^ 

H- ^—£1 • nv 2 ■ n£i • n 


— Si ■ n£i • 02 

8 £i • Fi£i • V 2 - Sfvl - {Si ■ V 2 y - 32£i • nSi ■ viV 2 ■ n 

2 miG L 

— 32£i • fivi ■ nSi ■ V 2 + 8 £i • nv 2 ■ nSi ■ V 2 — I 6 F 1 • F 2 (£i • n)^ 

+ 4?;|(£i • ay + 96Fi • nv 2 ■ n(£i • n)^ — 12(£i • n)^(F 2 • n)^ 

4G^C'iES'2m2 


(3.43) 


+ 


miT'^ 

G‘^Gies^'^2 

4^4 


£1 • n£i • F 2 + £1 • nSi ■ V 2 — 4£i • nSi ■ nv 2 ■ n 


(3.44) 


sivi - 12£i • Fi£i • V2 - 4£fFi • V2 + 2(£i • Fi)^ 


— 12£i • nvi ■ nSi ■ vi + 16£^Fi • nv 2 ■ n + 48£i • nSi ■ F 1 F 2 • n 
+ 48£i • nvi ■ nSi ■ F 2 — 5v‘f{Si ■ ay + 24Fi • F 2 (£i • n)^ — £1 (Fi • ay 
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(al) 


(a2) 


(bl) 



Figure 4. NNLO spin-squared cubic in G Feynman diagrams with no loops. 

G'^CiEg2m2 


— Iddui • nv 2 ■ n(5i • nf + 15(5i • nfivi ■ fty 


Si ■ nSi ■ vi 


-|- Si ■ nSi ■ vi — 2Si ■ S 1 V 2 ■ n — 35i • nSi ■ V 2 — • nSi ■ V 2 

— 2Si ■ nSi ■ nvi ■ n + 12Si ■ nSi ■ nv 2 ■ n (3.45) 

3.3 Cubic in G interaction 
3.3.1 Three-graviton exchange 

There are 4 diagrams at order with three-graviton exchange in this sector, which can be 
seen in figure 4. They do not contain linear in spin couplings, but rather up to three-graviton 
spin-squared ones. 

These diagrams are evaluated by: 


Fig. 4(al) 


C'ie52 


2 r 

Fig. 4(a2) =4 C'i£;s2 


- 3 5i • re 


G^rrel 


Sl-^[Si-n 


Fig. 4(bl) = 


Fig. 4(b2) = 


Gies"^ G^mim2 
4 r® 


Sf - 3 (Si-n 


St - 3 (Si-n 


(3.46) 

(3.47) 

(3.48) 

(3.49) 


4 rreir 

3.3.2 Cubic self-interaction with two-graviton exchange 

There are 4 one-loop diagrams at order G^ in this sector, which can be seen in figure 5. 
Two-graviton spin-squared couplings also appear here. 

These diagrams are evaluated by: 


Fig. 5(al) = - ^ rnim 2 


Fig. 5(a2) = - 8 


G^m^ 


Si ■ ft 


Fig. 5(bl) = — 


CiES'^ 


2 rreir' 
2 


5? + 5i • re 


Fig. 5(b2) =Gies2 ^ 




55f - 3 5i • re 


(3.50) 

(3.51) 

(3.52) 

(3.53) 
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Figure 5. NNLO spin-squared cubic in G one-loop Feynman diagrams. 



Figure 6. NNLO spin-squared two-loop Feynman diagrams. 


3.3.3 Two-loop interaction 

There are 11 two-loop diagrams at order in this sector, which can be seen in figure 6. 
Here, the cubic vertices have no time dependence, as in the spinl-spin2 sector [21], but unlike 
the spin-orbit sector [22]. This is similar to what happens in one-loop diagrams at the NLO 
level [16, 30, 33]. 

The irreducible two-loop tensor integrals, which are required here up to order 4, are 
reduced as in Appendix A of [21], see eqs. (A12), (A13), and in Appendix A of [22], see 
eq. (A2). 

These diagrams are evaluated by: 


Fig. 6(al) =0, 

Fig. 6(a2) =0, 

u- ^ G^mim2 

Fig. 6(bl) =—-.- 


35*^ -|- 5 I S'! • n 


Fig. 6(b2) = — 


1 G^mim2 
70 ^ 


135f - 25 5i • n 


(3.54) 

(3.55) 

(3.56) 

(3.57) 
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Fig. 6(b3) 
Fig. 6(b4) 
Fig. 6(b5) 


CiES'^ 


2 mir 
3C'ie52 G^mim2 

7 

Cies 2 G^mim2 


3Sf - 5 5i • n 


- 3 Si • n 


- 3 Si • n 


Fig. 6(cl) =0, 

Fig. 6(c2) =G,es^^ 



15 ( Si• n 


Fig. 6(d) =0, 


Fig. 6(e) 


4 G^mim2 
^ r5 


2Sl + 15 {Si • n) 



(3.58) 

(3.59) 

(3.60) 

(3.61) 

(3.62) 

(3.63) 

(3.64) 


Notice that the diagram in figure 6(cl) is equivalent to 0, despite its being irreducible, rather 
than of the factorizable kind. 


4 Next-to-next-to-leading order spin-squared potential 


The NNLO spin-squared potential for a binary with compact spinning components is split 
according to the number and order of higher-order time derivatives: 

.. (0) (1) (2) (3) 

^NLO + V + V -\- V ■ (4.1) 


The first part is independent of higher-order time derivatives and is thus an ordinary poten¬ 
tial, 


( 0 ) 

V= - 


IQrriir^ 


8Si • FiUiSi • V2 


+ 12Si • iJiSi • V 2 V 2 — 22S‘lvi ■ V 2 V 2 


42SiUii;i • V2 — 24Si • uiSi • V 2 V 1 ■ V 2 + llSiUiuf 
12^1(51 • Fi)^ 24Fi • F2 (Si • Fi)^ - 10r;2(Si • Fi)^ 


-I- 2ui(Si • V2)^ + 26Si (Fi • V2)^ + 19SiUi + + 42Si • FFi • nSi ■ vivf 


— ISSfvi ■ nv\v2 ■ n — 24Si • FSi • viv\v2 ■ n — 36Si • nvi ■ nv\Si ■ F2 

— 24Fi • FSi • F1F2 • FSi • V2 — 12Si • nv\v2 ■ nSi ■ F2 — 72Si • FFi • nSi • FiFi • F2 


— I2S1F1 • nv2 ■ nvi ■ V2 + 72Si • nSi • F1F2 • FFi • F2 -|- 72Si • FFi • FSi • F2F1 • F2 


-I- 24Si • FFi • nSi ■ Fiu^ -I- ISSfFi • nv2 ■ nv2 — 36Si • nSi • F1F2 • nv2 — 36Si • nFi • nSi • V 2 V 2 
+ 30uiFi • F2 (Si • nf — 9r;iu|(Si • F)^ -|- 30Fi • V2V2{Si ■ nf — 6(Fi • F)^(Si • Fi)^ 

— 36Si?;i(Fi • n)^ -|- 24Fi • nv2 ■ F(Si • Fi)^ -|- 12Si • FiSi • F2(Fi • n)^ -|- 6OS1F1 • F2(Fi • F)^ 

— 27Siu|(Fi • F)^ - 6(Si • Fi)^(F 2 • F)^ -h 9Sir;i(F2 • F)^ - 6(Fi • F)^(Si • V2)‘^ 

— 6(Si • F)^(Fi • F2)^ — 21(Si • n)^vf — 21(Si • n)‘^V 2 + 90Fi • Fr;iF2 • F(Si • F)^ 

— I8OF1 • FF2 • FFi • F2 (Si • F)^ -|- 90Fi • FF2 • nv^iSi ■ F)^ -|- 6OS1F2 • F(Fi • F)^ 

— 6OS1 • FSi • F1F2 • F(Fi • F)^ -|- 6OS1 • nv2 ■ FSi • F2(Fi • F)^ -|-15^1 (Si • F)^(Fi • F)^ 

— 45Si (Fi • F)^(F2 • F)^ -b 6OS1 • FFi • FSi • Fi(F2 • F)^ -b 15r;f (Si • n)‘^{v2 ■ F)^ 
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mir^ 
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+ 445i • nv2 ■ nSi ■ V2 — Vi{Si ■ n)^ + 21;i • V2{Si ■ n)^ + 21;!(5i • n)^ — 365i(l72 • n)^ 
+ 241;i • nv2 ■ n{Si ■ n)^ + 6(5i • n)^(l7i • n)^ — 12(5i • n)‘^{v2 ■ 77)^ 


G‘^m2 

47-4 
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2Slvi ■ V2 


35^2 
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G'^CiES^fn2 
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G^mim2 ' 

147-5 
21G^ (71^52 7712 ■■ 

2r5 


Sl + 25(5i • nf 


+ 
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(5i • nf - 


Sj - 3(5i • nf 


IdG^CiEs^nT-^ 

4mir® 


23G^C'i£;52mim2 

2^ 


Sf - 3(5i • nf 
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(4.2) 


The part linear in higher-order time derivatives is further split due to its length as 

( 1 ) 


V=Va + Vs, 


(4.3) 


where 
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SSfvfSi ■ n — 16vi ■ nSi ■ viSi ■ Si + 165il7i • nvi ■ Si + 28SiVi ■ S 1 V 2 • n 
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mir^ 

G‘^Cies^i^2 


Si ■ nSi ■ Si — 115fa2 • n + 145i • nSi ■ 02 — Si ■ n{Si ■ Sy 


SfSi ■ n + 2Si ■ n{Si ■ Sy 


G‘^m2 


SSfSi ■ n — 3Si ■ nSi ■ Si + SIS 2 ■ n — Si ■ nSi ■ S 2 + (1 -H- 2), 


(4.4) 


and 


V- = 


GGieS‘^1^2 

4mir2 


4ni • nSi ■ viSi ■ vi — 4S'i • • nvi — 2Si ■ viSi ■ V 1 V 2 ■ n 


— 3Si ■ Sivfv 2 ■ n — 5ni • nSi ■ viSi ■ V 2 + Si ■ • V 2 + 2Si ■ V 1 V 2 ■ nSi ■ V 2 

— 5vi ■ nSi ■ viSi ■ V2 + Si ■ nvfSi ■ V 2 + 2Si ■ V 1 V 2 ■ nSi ■ V 2 + 6ni • nSi ■ V 2 S 1 ■ V 2 
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35i • SiVi ■ n — GSi • SSi ■ vi — 5Si ■ nSi ■ vi + 2Si ■ S 1 V 2 ■ n + 2Si ■ nSi ■ V 2 


+ 2Si ■ nSi ■ V 2 + 145i • n5i • nvi ■ n — 65i • nSi ■ nv 2 ■ S 
G‘^Gies^'^2 


+ 


2r3 
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— ^ „ —- 2Si ■ SiVi ■ n — 3Si ■ nSi ■ vi — 3Si ■ nSi ■ vi + 305i • S 1 V 2 ■ S 

2miG L 
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+ (lo2). 

The contribution with two higher-order time derivatives reads 


(4.5) 
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+ 135i • Si{v 2 ■ n)^ — 4^1 • nv 2 ■ nSf — 3Si ■ nSi ■ n{v 2 ■ n)^) + (25i • ai5i • 02 + 135fai • 0,2 
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— Si ■ fiSi ■ vid2 ■ fi — Si ■ nSi ■ vid2 ■ fi + 35i • fivi ■ fiSi ■ 3,2 + 35i • fivi ■ fiSi ■ 3,2 


+ 7Si ■ nSi ■ nvi ■ 02 — 65i • nSi ■ nv 2 ■ 02 + 2ai • nv 2 ■ nSf + 35i • n5i • nvi ■ no 2 • n) 

— 2[8Sivl + 8^1 • • V 2 — 8Slvi ■ V 2 + 29S‘fv2 — 8{Si ■ vi)"^ + 2{Si ■ V 2 )^ + 8 Sfvi ■ nv 2 ■ n 

— 8Si ■ nSi ■ V 1 V 2 ■ n + 125i • nv 2 ■ nSi ■ V 2 — “ 13Sf{v2 ■ n)^ 


+ 3{Si ■ n)‘^{v2 ■ n)^) 
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G'^CiEs^rn^ " 
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2mir^ 

G'^Cies'^'^2 


{Si -81 + Si- nSi • n) + {Sf + {Si • n)^) 
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(4.6) 


Finally, several terms with three higher-order time derivatives appear in the potential as well, 

4 i;2 • nSf (l3Si • 5i02 • re — 7Si ■ nSi ■ 02 — 7Si ■ nSi ■ 02 


(^)_ ^GCiES2m2 


8mi 


— Si ■ nSi ■ 332 ■ 3 '^ + (135^02 • re — 14,Si • 3Si • 02 — ^2 • 3{Si ■ re)^) -|- (1 -H- 2). (4:.7) 


As explained in [16], since this potential contains higher order time derivatives of the variables, 
the EOM, which are derived from it via a variation of the action, also contain higher order 
time derivatives, which must be eliminated, using lower order EOM. In a forthcoming work we 
are going to reduce these higher order time derivatives at the level of the potential, following 
the procedure of redefinition of variables, outlined in [29, 37], where it was extended for spin 
variables in [29]. We are going to provide the EOM, and to transform this potential to a 
Hamiltonian, from which gauge-invariant relations for the binding energy can be obtained in 
a straightforward way. These relations are of great value for refining, for instance, the EOB 
Hamiltonian, and the gravitational waveforms derived from it. 


5 Conclusions 

In this work we derived for the first time the NNLO spin-squared interaction potential for 
generic compact binaries via the EFT for gravitating spinning objects in the PN scheme, 
which was formulated in [16]. This correction, which enters at the 4PN order for rapidly 
rotating compact objects, completes the conservative sector up to the 4PN accuracy. This 
high level of PN accuracy is necessary in order to detect gravitational waves. The robust¬ 
ness of the EFT for gravitating spinning objects is shown here once again, as was already 
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demonstrated in [16, 18, 22], which obtained all spin dependent sectors, required up to the 
4PN accuracy, and was initiated at the NNLO level in [21, 29]. The EFT of spinning objects 
allows to directly obtain the EOM from the potential via a proper variation of the action 
[16, 29], and to obtain corresponding Hamiltonians in a straightforward manner [16]. In a 
forthcoming paper we are going to derive these for the NNLO spin-squared potential obtained 
in this paper. 

The spin-squared sector is an intricate one, as it requires the consideration of the point 
particle action beyond minimal coupling, and involves mainly contributions with the spin- 
squared worldline couplings, which are quite complex, compared to the worldline couplings, 
that follow from the minimal coupling part of the action. Yet, there are also contributions 
to this sector, involving the linear in spin couplings, as we go up in the nonlinearity of the 
interaction, and in the loop order. Thus, there is an increase in the number of Feynman 
diagrams and in complexity with respect to the NLO, which is even more excessive than 
that, which occurs in the spinl-spin2 sector [16, 21, 30]. We have here 64 diagrams, of which 
more are higher loop ones, where we have 11 two-loops, and tensor integrals up to order 6. 
The application of the “NRG” helds [31, 32], and the gauge choices of the EFT for gravitating 
spinning objects [16] then help to render the computations more efficient. Lastly, we checked 
most of the computations, using Mathematica. 

Finally, it should also be noted, that the relation between the Wilson coefficients, ap¬ 
pearing in this sector, and in [18], required for generic compact objects at the 4PN accuracy, 
and the multipole moments used in numerical codes, e.g. [38-40], as well as the relation with 
the Geroch-Hansen multipoles for black holes [41, 42], should be worked out via a formal 
EFT matching procedure. Glearly, this involves subtleties [39], and is left for fntnre work. 
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